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ABSTRACT 

Context. Some recent observations of the delay-time distribution (DTD) of Type la supernovae (SNe la) seem to uphold the double degenerate 
^0 ° (DD) scenario as the progenitor model of SNe la, but the core-degenerate (CD) scenario remains a strong competitor to the DD one. 
(-H . Aims. We investigate the effects of metallicity and the different treatments of common envelope (CE) on the DTD of SNe la by considering the 
f») ' DD and CD scenarios, and check the suggestion that the total mass of DD system is the main dependent variable of Phillips relation, 
l Methods. We perform a series of Monte Carlo simulations based on a rapid binary evolution code and consider two treatments of CE evolution, 
O ■ i.e.Q'-formalism and y-algorifhm. 

■j-J ' Results. We find that only when the ff-formalism is considered with a high CE ejection efficiency, may the shape of the DTD for DD systems 
^ , be consistent with that derived observationally, i.e. a power law of ~ t -1 , while the value of the birth rate of SNe la marginally matches 

i i observations. For the fy-formalism with a low CE ejection efficiency and the y-algorifhm, neither the shape of the DTD nor the value of the 

birth rate can be compared with those of the observations. Metallicity may not have a significant influence on the shape of DTD, but a lower 
metallicity may lead to a slightly higher birth rate of SNe la by a factor of 2, especially for SNe la with long delay times. If the results for the 
single degenerate (SD) channel are incorporated into those for the DTD, both the shape of DTD and its value may be closely consistent with 
observations for SNe la younger than 2.5 Gyr, and SD and DD channels provide comparable contributions to the total SNe la, while for SNe la 
with delay times longer than 2.5 Gyr, DD is the dominant channel and the birth rate is lower than that derived from observations by a factor up 
to about four. In addition, we calculate the evolutions of various integral parameters of DD systems, and do not find any one suitable to explain 
the correlation between the brightness of SNe la and its delay time. Moreover, there are three channels producing core-degenerate(CD) systems 
in that may contribute a few SNe la, but the contribution of CD systems to the total SNe la is no more than 1%. 
' Conclusions. There may be other channels or mechanisms contributing to SNe la with long delay times. 
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Vh , 1. Introduction Leibundgut 2000; Pailhasarathy et al. I20071 Wang & Han 

=3 |2009b). There is a consensus that SNe la result from the ther- 

Type la supernovae (SNe la) have been found to be important , , • e u rr^rw u-, a 

Jr F v ' F monuclear explosion of a carbon-oxygen (CO) white dwarf 

to many astrophysical fields, especially as principal distance in- (WD) fa a binafy system (Hoyle & Fowler[rm According to 

dicators to measure cosmological parameters, where their use , f . u f ., .. , lrrv , 

fe F the nature of the companions of the mass-accreting WDs, two 
resulted in the discovery of the accelerating expansion of the U1 .. e CXT T , , ■• , 

i ^ f possible progenitors of SNe la have been discussed over the 

Universe (Riess et al. 1998; Schmidt et al. 1998; Perlmutter et . ■ ■ • • , , , /C t->\ a ^ u 
v 1 ' ' ' past three decades, i.e. the single degenerate (SD) model where 

al. 1 19991 >. This result was extremely exciting and implied that • • , . ,,,,, , „ TU i m-ni vr 

1 ' the companion is a normal star (Whelan & Iben 1973 , Nomoto 

dark energy exists. At present, SNe la are regarded as critical et ^ rjgg^ and the double degenerate (DD) mode i where a 

cosmological probes for testing both the evolution of the dark CQ WD merges ^ anQther CQ WD (Iben & Tutukoy 19g4; 

energy equation of state with time and the evolutionary his- , , ■ , i lr , 0/ |i A . . , . , ij t 

bJ M £ Webbink |1984| l. At present, both models cannot be ruled out 

tory of the universe (Riess et al. [2007] Kuznetsova et al.lffll completely by observations (see the review by HowellEHID. 

Howell et al. 1200951 F 3 3 J 

However, the precise nature of SNe la remains unclear, To distinguish between these different progenitor models, 
especially their progenitor (Hillebrandt & Niemeyer EOQO] it is important to measuring the delay-time distribution (DTD), 
where the delay time is the elapsed time between the primor- 
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event and the DTD is the SN birth rate versus their delay time 
for a single starburst. An increasing amount of observational 
studies of the DTD of SNe la have found that the DTD fol- 
lows a power-law form of r l (Totani et al. [2008, Maoz & 
Mannucci l201 II) . which is difficult to interpret in terms of the 
SD model (Meng & Yang 2010a), but in contrast can be natu- 
rally explained by the DD model (Yungelson & Livio 2000; 
Mennekens et al. 2010). Some theoretical and observational 
studies have shown that metallicity may have a strong influ- 
ence on the birth rate of SNe la (Khan et al. 1201 11 Meng et 
al. 1201 lab . We wish to establish the effect of metallicity on the 
birth rate of SNe la in the DD scenario. 

When SNe la are used as distance indicators, the Phillips 
relation is adopted, which is a linear relation between the ab- 
solute magnitude of SNe la at maximum light and the mag- 
nitude drop in the B band light curve during the first 15 days 
following the maximum light (Phillips 1993). This relation was 
motivated by the observations of two peculiar SN events (SN 
1991bg and SN 1991T) and implies that the brightness of SNe 
la is determined mainly by one parameter. The amount of 56 Ni 
formed during the SN explosion dominates the maximum lu- 
minosity of SNe la (Arnett 1982). However, the origin of the 
variation in the amount of 56 Ni for different SNe la is still un- 
clear (Podsiadlowski et al. 12008'). Observationally, the most lu- 
minous SNe la always occur in spiral galaxies, while both spi- 
ral and elliptical galaxies host systematically dimmer SNe la, 
which leads to a dimmer mean peak SN brightness in elliptical 
than spiral galaxies (Hamuy et al. 1996; Brandt et al. [2010). 
In addition, the mean peak brightness of SNe la in a galaxy 
varies less in the outer than the inner regions (Wang et al. |1997l 
Riess et al. |1999l ). In other words, age could be the most impor- 
tant factor in determining the luminosity of SNe la, and dim- 
mer SNe la have a wide age distribution (Gallagher et al. 2008 ; 
Neill 2009; Howell et al. 2009b). From these observations, we 
concluded that the range and average value of the brightness 
of SNe la decrease with their delay time. If the maximum lu- 
minosity of SNe la is mainly determined by one parameter as 
shown by the Phillips relation, we may expect that the range of 
the parameter decreases, and its average value either increases 
or decreases with the age of SNe la. 

Many efforts have been made to resolve this problem. 
Some multi-dimensional numerical simulations have shown 
that varying the ignition intensity at the center of WDs or 
the transition density from deflagration to detonation might 
enable us to determine the underlying cause of the observed 
correlation between peak luminosity and light-curve width 
(Hillebrandt & Niemeyer 120001 Hoflich et al. 120061 I20T01 
Kasen et al. 120101 ). In addition, the ratio of nuclear-statistical- 
equilibrium to intermediate-mass elements in the explosion 
ejecta may be a key parameter in determining the width of SN 
la light curve and its peak luminosity (Pinto & Eastman 2001 1 
Mazzali et al. l2lMl 120071) . Lesaffre et al. d2006l ) suggested that 
the central density of the WD at ignition may be the origin 
of the Phillips relation in a systematic study of the sensitivity 
of ignition conditions for H-rich Chandrasekhar-mass single- 
degenerate exploders to various properties of the progenitors 
(see also Podsiadlowski et al. 120081 ). which was upheld by de- 
tailed multi-dimensional numerical simulations of explosion 



(Krueger et al. 2010). Moreover, metallicity may also affect the 
production of 56 Ni, thus the maximum luminosity, both in the- 
ory (Timmes et al. 2003, Travaglio et al. 120051 Podsiadlowski 
et al. 120061 Bravo et al. 120 1 Oal ) and observations (Branch & 
Bergh 119931 Hamuy et al. fl996l Wang et al. 119971 Cappellaro 
et al. 119971 Shanks et al. 120021 ). On the basis of the suggestion 
that the ratio of the carbon to oxygen (C/O) of a white dwarf at 
the moment of explosion could be the dominant parameter de- 
termining the production of 56 Ni (Nomoto et al. 119991 120031 ). 
Meng & Yang d201 II) found that the effects of metallicity and 
C/O on the production of 56 Ni may complement each other. 
Most of these discussions above focused on the Chandrasekhar 
mass model, in which the WDs explode as SNe la when their 
masses are close to the Chandrasekhar mass limit. However, 
the sub-Chandrasekhar mass model may still be able to explain 
SNe la such as 1991bg, as well as normal SNe la (Sim et al. 
120101 Ruiter et al. 1201 II ). In addition, observations of several 
very bright SNe la implied that their progenitor WDs might 
have a super-Chandrasekhar mass (Astier et al. 120061 Howell 
et al. 120061 Hicken et al. 120071 Scalzo et al. l20T01 Yuan et al. 
[20101 Tanaka et al. I20T01 Yamanaka et al. I20T01 ). Although 
these very bright SNe la may also be explained in terms of 
the SD scenario (Chen & Li l2009l Liu et al. l2009l Hachisu et 
al. 12012), that the DD systems are the progenitor of the very 
bright SNe la still cannot be completely ruled out at present. 
Howell d201 II ) suggested that the DD scenario has a natural ex- 
planation of the higher SN luminosity in young environments, 

1. e. younger, more massive stars produce more massive white 
dwarfs that have more potential fuel than less massive mergers, 
and then produce a bright SNe la (see also Maoz & Mannucci 
1201 II ). This suggestion still needs to be checked carefully. 

However, earlier numerical simulations showed that the 
most probable fate of the DD merger is an accretion-induced 
collapse (AIC) and, finally, neutron star formation (see the re- 
view by Hillebrandt & Niemeyer 12000) . Even if the merger 
survived the AIC, a super wind from giant-like structure before 
supernova explosion would occur and the remnant might lose 
about 0.5 M Q and shrink in mass to below the critical mass for 
explosion (Soker 1201 11 . Following the suggestion of Sparks 
& Stecher JT9741 ) and Livio & Riess d2003l . Kashi & Soker 
d201 II) developed a core-degenerate (CD) model to overcome 
the drawbacks of the DD model, in which a white dwarf merges 
with the core of an asymptotic giant branch (AGB) star shortly 
after a common envelope (CE) phase. However, they did not 
consider the evolution of the birth rate from the CD supernova. 
In this paper, we study the evolution in terms of a detailed bi- 
nary population synthesis (BPS) study. 

In Sect. [2] we simply describe our method, and present the 
calculation results in Sects. |3]and|4] In Sect. [5] we show our 
discussions and main conclusions. 

2. METHOD 

2.1. Common envelope 

The CE phase is very important to the formation of DD 
systems. During binary evolution, the mass ratio (q - 
Mtonoi-M^accretor) is a crucial parameter. If it is larger than a crit- 
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ical mass ratio, q c , mass transfer between the two components 
is dynamically unstable and a CE forms (Paczynski 1976 ). The 
ratio q c varies with the evolutionary state of the donor star 
at the onset of Roche lobe overflow (RLOF) (Hjellming & 
Webbinl JT987l Webbink [19881 Han et al. 120021 Podsiadlowski 
et al. [20021 Chen & Han 120081 . In this study, we adopted q c 
= 4.0 when the donor star is on the main sequence (MS) or 
crossing the Hertzsprung gap (HG). This value is supported by 
detailed binary evolution studies (Han et al. 2000; Chen & Han 
2002 2003 ). If the donor star is on either the first giant branch 
(FGB) or the AGB, we use 



1c 



1.67 



\ M . 



/2.13, 



(1) 



where M c is the core mass of the donor star, and x = 
dln/?i/dlnM is the mass-radius exponent of the donor star, 
which varies with composition. If the mass donors are naked 
helium giants, then q c = 0.748 based on Eq. ([TJ (see Hurley et 
al.|2002]for details). 

Embedded in the CE are the dense core of the donor star 
and the secondary. Owing to frictional drag with the envelope, 
the orbit of the embedded binary decays, and a large part of the 
orbital energy released in the spiral-in process is injected into 
the envelope (Livio & Soker ll988T ). 

The CE evolution is very complicated and different authors 
may use different methods to treat it in their BPS studies. There 
are two dominant methods for treating the CE evolution, i.e. 
the of-formalism ensuring energy balance and the y-algorithm 
ensuring angular momentum balance. The a-formalism, which 
is widely used, may closely reproduce the orbital-period dis- 
tribution of WD + MS systems as noted by Zorotovic et al. 
(120101 ) (see also Hurley et al.[2002]and Webbink|2007]), while 
it may be unable to describe the production of a close pair of 
white dwarfs. To solve this problem, Nelemans et al. ((2000) 
and Nelemans & Tout ( 2005 ) developed the y-algorithm, which 
may explain the formation of all kinds of close binaries. In this 
paper, both of these descriptions are applied. 

For the cr-formalism, the final orbital separation cif after CE 
phase is obtained by 



af M c / + 2M e flj 



a, 



acE^mR 



(2) 



where a; is the initial orbital separation at the onset of the CE, 
the masses M, M c , and M e are those of the donor, the donor 
core, and the envelope , respectively, Ri is the radius of the 
donor, and m is the companion mass. The parameter cce is the 
CE ejection efficiency, i.e. the fraction of the released orbital 
energy used to eject the CE, and A is a structure parameter re- 
lying on the evolutionary stage of the donor. At present, the 
value of acE is very uncertain and may vary from 0.4 to 3.0 
(see the review of Ivanova l201 II ). Because the thermal energy 
in the envelope is not incorporated into the binding energy in 
this paper, acE mav t> e greater than 1 (see Han et al. 1 1995 1 for 
details about the thermal energy). In this paper, we set cuce to 
either 1.0 or 3.0 when studying the DD and CD scenarios, and 
acE = 1.5 or 2.0 were also tested. For A, we assumed that it is 
a constant (A - 0.5, de Kool, van den Heuvel & Pvlyser fT987l ). 



In the y-algorithm, Of was obtained from 



a 1= (M\ (Mo±m\L_ J^_\ 
fli \M C \M + ml\ ■ M + ml 



(3) 



where jce is a free parameter. Based on the results of Nelemans 
& Tout (120051 ). we set y CE = 1.5. 

2.2. Basic parameters of our Monte Carlo simulations 

To investigate the birthrate of SNe la in DD systems, we fol- 
lowed the evolution of 1 x 10 7 binaries via Hurley's rapid bi- 
nary evolution code (Hurley et al. 2000 2002 ). The primordial 
binary samples are generated in a Monte Carlo way and a cir- 
cular orbit was assumed for all binaries. The basic parameters 
for the simulations were as follows: (1) a constant star forma- 
tion rate (SFR) of 5 M Q yr over the past 15 Gyr or a sin- 
gle starburst of 10 n M Q ; (2) the initial mass function (IMF) of 
Miller & Scalo (119791 ): (3) a constant mass-ratio distribution; 
(4)all stars are members of binary systems and the distribution 
of separations is constant in log a for wide binaries, where a 
is the orbital separation, and falls off smoothly at small sepa- 
rations, where a = 10/? Q is the boundary for wide and close 
binaries; (5) metallicities were chosen to be Z = 0.02, 0.002, 
and 0.0001 (see Meng & Yang 1201 Oil for details of these input 
parameters). 

We tracked the evolutions of the 1 x 10 7 binaries until they 
become DD systems. Following the DD system, gravitational 
wave radiation (GWR) dominates the evolution of the system 
on a timescale few (Landau & Lifshitz [19621 ) 



?Gw(yr) = 8 x 10 7 x 



(M l+ M 2 )^ rm 
M\M 2 



(4) 



where P is the orbital period of the DD system in hours, and 
Mi and M 2 are the masses of the two WDs in solar masses, re- 
spectively. The time elapsed from the birth of primordial binary 
system to the occurrence of SN la was then assumed to be equal 
to the sum of the timescale on which the primordial secondary 
star becomes a WD and the orbital decay time. We assume that 
if the total mass of a DD system is M t = Mi + M 2 > 1.378M Q 
and the elapsed time is shoter than 15 Gyr, a SN la is produced. 

2.3. Evolution channel for the CD system 

A system consisting of a WD and an AGB core survives a CE 
phase and then merges shortly after. The core is more massive 
than the WD and when the merging process starts, it is still very 
hot and has a larger radius than the cold WD. Hence, the core 
may be destroyed to form an accretion disk around the WD. 
Since the less massive WD has a shallower gravitational po- 
tential, the temperature in the accretion disk does not reach the 
ignition temperature and AIC can be avoided. This scenario is 
only applicable when the merging process occurs within ~ 10 5 
yr of the CE phase. The merger is unlikely to explode soon un- 
til a large part of the angular momentum is lost by means of the 
magneto-dipole radiation torque (see Ilkov & Soker l20T2l for 
details). Thus, as mentioned in Soker (201 1), three key factors 
distinguishing the CD scenario from the DD channel are that 
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in the former scenario: (1) the hot core is more massive than 
its cold companion WD, (2) the CD system will merge within 
~ 10 5 yr after the CD system formsQ, and (3) its delay time is 
mainly due to the spinning-down time of the merger product, 
which results from a magneto-dipole radiation torque (Ilkov & 
SokerlUn). Ilkov & Soker (I20T21 and Soker (EOTTt only in- 
troduced one channel to produce CD systems. Three channels 
could actually produce CD systems, fulfilling the above crite- 
ria. However, different treatments of CE may trigger different 
channels as follows. In the following description, the primary 
is the more massive star in a primordial binary system, while 
the secondary is the less massive one. 

Case 1: (wind + RL + CE) This channel can be encoun- 
tered by both a-formalism and y-algorithm, but is more com- 
mon for y-algorithm, and is similar to that described in Soker 
( 1201 11 1 and Ilkov & Soker (2012), which includes two sub- 
channels. For sub-channel A, the primordial zero-age main se- 
quence (ZAMS) mass of a primary is in the range from 4.0 M 
to 5.0 M Q , and the mass ratio {m%jm\) is relatively high(0.85 - 
0.95). The primordial binary system has a very wide separation 
(wider than 2300 R Q ), which permits the primary to evolve into 
the thermal pulsing asymptotic giant branch (TPAGB). Before 
the primary fills its Roche lobe, it loses a lot of material by 
its wind, which results in the following stable RLOlU At this 
stage, the secondary is still a MS star. The system becomes a 
CO WD + MS system after the RLOF. The MS secondary then 
has a mass higher than the primordial mass of the primary. The 
WD + MS system continues to evolve and the MS secondary 
fills its Roche lobe when it becomes a AGB star. Because of the 
large mass ratio at this stage, the system enters a CE phase. If 
the CE can be ejected, a CD system forms and merges within 
~ 10 6 yr. 

For sub-channel B, the primordial primary usually has a 
mass of 3.0 - 3.5M , and the mass ratio is higher than 0.99. 
However, the primordial separation for sub-channel B is not 
as wide as that for the sub-channel A (18O(M? - 21O(M? ). 
The evolutionary path of the sub-channel B is much similar 
to that of the sub-channel A and the difference is that when the 
primary becomes an AGB star, the secondary is a horizontal 
branch (HB) star (central helium burning), not a MS star. 

Case 2: (RL + 4CE) This channel is only encountered when 
y-algorithm is adopted. The primordial primary has a mass 
higher than 6.5 M and the mass ratio is larger than 0.94, even 
close to 1, but the primordial separation is only about 100 R Q . 
Owing to the small separation, when the primary crosses the 
HG, it fills its Roche lobe and a stable RLOF occurs, which re- 
sults in the secondary being more massive than the primordial 
primary. The separation of the components then increases dur- 
ing the stable mass transfer, the primary continues to evolve to 
the FGB, and the mass transfer then becomes unstable, leading 
to the first CE. After the CE ejection, the system consists of a 



1 Note that the time here is the GWR timescale after the CD system 
forms and is not the delay time of a SN la. 

2 The treatment of RLOF in the Hurley's rapid binary evolution 
code is a substantially revised version of that presented by Tout et 
al. (Q997 ) and the radius-mass exponent £ defined by Webbink ( 1985 ) 
is used to describe the stability of mass transfer. 



naked helium star and a MS star. According to the y-algorithm, 
the separation does not shrink greatly after the first CE ejection. 
Shortly afterwards, the MS star evolves to the FGB and fills its 
Roche lobe where the mass transfer is dynamically unstable, 
and the system then enters into the second CE phase. After the 
CE ejection, the secondary also becomes a naked helium star, 
and the system consists of two helium stars. The mass of the he- 
lium star from the secondary is higher than that of the primary 
because the secondary at the MS stage is more massive than the 
primordial primary for prior stable RLOF (mi e [1.22, 1.35]M 
and m.2 e [1.4, 1.6]M , respectively). The helium star from the 
primary firstly exhausts its central helium and then begins to 
cross the helium HG. If it fills its Roche lobe, the third CE is 
expected for a large mass ratio. However, the mass of the third 
CE is low compared to the two prior CEs, and the third CE 
evolution cannot shrink the separation significantly. The CE is 
again ejected and then the primary becomes a CO WD. Helium 
is also finally exhausted at the center of the second helium star. 
As the second helium star crosses the helium HG, it also fills 
its Roche lobe, leading to the fourth CE. After the ejection of 
the last CE, a CD system forms. Because of the four CE evo- 
lutions, the separation of the CD system is so close that the 
system merges within ~ 10 6 yr. 

Case 3: (RL + CE + RL/wind + CE) This channel is en- 
countered when an a-formalism is adopted. The primordial pri- 
mary has a mass of 4.3 -5.2M and its mass ratio is lower than 
in the above two cases (0.68 - 0.78). The range of primordial 
separations is from 20 R Q to 55 R Q . Because of the small sep- 
aration, the primary fills its Roche lobe as it crosses the HG, 
where a stable mass transfer occurs. The primary loses its en- 
velope and becomes a naked helium star. At the same time, the 
secondary accretes a lot of material, becoming more massive 
than primordial primary. Owing to the stable RLOF, the sep- 
aration greatly increases, and the secondary may then evolve 
to its HG/RGB phase and fill its Roche lobe. The subsequent 
mass transfer is dynamically unstable and a CE forms. After 
the CE ejection, the system consists of two naked helium stars, 
and owing to the prior stable RLOF, the helium star evolving 
from the secondary is more massive (m\ e [O.7,O.85]M and 
m.2 6 [1.25, 1.6]M Q , respectively). The helium star that evolves 
from the primary firstly exhausts its central helium and begins 
to cross the helium HG. At this stage, if the helium star fills its 
Roche lobe, a stable RLOF occurs and the helium star trans- 
fers its envelope onto the second helium star and then becomes 
a CO WD. It otherwise loses its envelope in terms of a wind 
to become a WD. The second helium star also finally exhausts 
its central helium, leading to the expansion of its envelope. Its 
Roche lobe is then filled and the system enters into the second 
CE phase for unstable mass transfer. After the CE ejection, a 
CD system forms and merges within 10 6 yr. 

In this paper, we assumed that (1) if a WD produced from a 
primordial secondary is more massive than its companion WD, 
the system is a potential CD system; (2) if the potential CD 
system merges after the emission of GWR within a fewxlO 5 
yr, the merger will explode as a SNe la. In addition, to study 
the birth rate of SNe la from CD systems, we traced the evolu- 
tion of 1 x 10 8 binaries. To check the influence of the merger 
timescale by GWR on the birth rate of SNe la in the CD sys- 
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log(«/yr) 

Fig. 1. Evolution of the birthrates of SNe la for a single star- 
burst of 1O U M , where Z = 0.02. Solid, dashed, and dot- 
dashed lines show the cases with uqe = 1-0, o;ce = 3.0, and 
7ce = 1.5, respectively. The dotted line represents a power-law 
function of form fit) oc r 1 . 

tem, we considered the cutting time as a free parameter (see 
SecLgJ. 

There is actually a "pollution" channel, that also fulfills our 
criterion (1) for the potential CD systems. In the channel, the 
more massive WD also evolves from secondary, but forms ear- 
lier than the second WD. The primordial primary has a mass of 
around 5M , and the mass ratio is close to 1. The primordial 
separation is in the range of 130 - 16O/? , which permits the 
primary to evolve to the FGB. A dynamically unstable mass 
transfer then results in a CE phase. Following the CE ejection, 
the system consists of a naked helium star and a MS one. On the 
basis of the y-algorithm for CE evolution, the orbital separation 
does not shrink greatly after the first CE evolution, which per- 
mits the secondary to evolve to the early AGB phase, and then 
the second CE. After the CE is ejected, the secondary becomes 
a WD and the separation shrinks greatly. After the helium star 
exhausts its central helium, the helium star also fills its Roche 
lobe, and the system enters into the third CE phase. After the 
CE ejection, the system consists of a less massive hot core and 
a more massive cold WD. The following GWR timescale for 
the channel may be as long as 1 .5 Gyr. 

3. Results for a normal DD channel 

3.1. Birth rate for a single starburst 

In Fig. Q] we show the evolution of the birthrates of SNe la 
from DD systems for a single starburst with Z = 0.02, where 
the evolution is characterized by a DTD. For comparison, we 
also show a power-law function of the form /(f) oc t by a 
dotted line in Fig. Q] We can see from Fig. Q] that the shape of 
the DTD depends on the treatment of CE. Some recent mea- 
surements measure a similar power-law shape with an index of 
-1 (Maoz & Mannucci l2011l) . In Fig.Q] only when a C E = 3.0 
(a'CE'i = 1-5), does the DTD follow a power-law shape of t . 
If cxce - 1.5, we may obtain a power-law DTD, while for 
acE = 1-0 (acE^ = 0.5), not only does the shape not follows 
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Fig. 2. Evolution of the birthrates of SNe la for a single star- 
burst of 1O H M with different metallicities, where acE - 3.0. 
The circle is from Mannucci et al. (2005), the filled circles 
are from Totani et al. (J2008 >, stars are from Maoz & Badenes 
(2010), squares are from Maoz et al. (2010 1, and triangles are 
from Maoz et al. (1201 11 1. 



a power-law shape, but the birth rate is also significantly lower 
than that for acE = 3.0 (see also Ruiter et al. I20091 >. This is be- 
cause that a low acE generally means a more heavy shrinkage 
of the orbital separation to eject CE. The system is then more 
likely to merge into a single star, not a DD system. However, 
the birthrate is almost unaffected by variations in the delay time 
in the y-algorithm, i.e. the DTD does not follow a power law. In 
the y-algorithm, the separation of a DD system surviving a CE 
phase is usually larger than that in the a-formalism (see Eq.0, 
which leads to a very long GWR timescale, and then a lower 
birth rate within 15 Gyr. It is noteworthy to point out that the 
y-algorithm favors the production of SNe la with a long delay 
time (Meng et al. 1201 lbV When oqe = 3.0, about two thirds 
of the SNe la explode within 1 Gyr, which is qualitatively con- 
sistent with observations (see the review of Maoz & Mannucci 
201 1 ), while in the other two cases, the SNe la older than 1 Gyr 
are more common. 

In Fig.Q] there is a single spike at early times, because the 
SNe la with delay times shorter than 10 8 yr are very rare for the 
DD channel. To save CPU time, we defined a binary sample to 
be 10 7 binaries. If we enlarge the sample to be large enough, the 
single spike should disappear and the line become continuous. 

We also check the dependence of the DTD shape on metal- 
licity, finding little significant influence using both the a- 
formalism and y-algorithm. Fig.|2]shows the effect of metallic - 
ity on the DTD of a C E = 3.0. We note that for the SNe la with 
the shortest delay time, a high metallicity leads to a systemati- 
cally earlier explosion time of SNe la, which is mainly due to 
the effect of metallicity on the stellar structure. For these SNe 
la, their delay time is dominated by the evolutionary time of the 
primordial secondary and the contribution of GWR is negligi- 
ble. For a certain evolutionary stage, a star with a high metal- 
licity generally has a larger radius, and may then fill its Roche 
lobe earlier. Hence, a given binary system with a high metallic- 
ity usually enters the CE phase earlier. If a DD system survives 
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Fig. 3. The evolution of the birth rates of SNe la for a constant 
star formation rate (Z=0.02, SFR=5M Q yr~ 1 ). Solid, dashed, 
and dot-dashed lines show the cases of c?ce = 1.0, acE = 3.0, 
andycE = 1.5, respectively. 









// «CE = 3 -° ; 




/ / Z=0.02 




I ■ Z-0.002 : 




'/ / Z=0.0001 ■ 




i 1 
i 

i i - 

1 i : 
, i , , i i , , , , i , , , , i , , , , i , , , , i , , , , ■ 



7.5 8 8.5 9 9.5 10 10.5 



log(i/yr) 

Fig. 4. The evolution of the birth rates of SNe la for a constant 
star formation rate (SFR=5M Q yr~ 1 ) with different metallicities, 
where acE = 3.0. Solid, dashed, and dot-dashed lines show the 
cases of Z = 0.02,0.002, and 0.0001, respectively. 



Table 1. DTD normalization(A^sN/A^«) from DD mergers for 
different metallicities and different treatments of CE in unit of 
10- 3 /M G . 



z 


ttCE =1.0 


Q-cE = 3.0 


7ce = 1.5 


0.02 


0.098 


0.58 


0.26 


0.002 


0.13 


0.87 


0.24 


0.0001 


0.22 


0.64 


0.29 



3.2. Birth rate for a constant star formation rate 

Fig. |3] shows the Galactic birth rates of SNe la (i.e. Z = 0.02 
and SFR= 5.OM yr~') from the DD channel. In the figure, 
the Galactic birth rate is around 0.5-3.0xl0 _3 yr , which is 
marginally consistent with that inferred from observations (3- 
7xl0~ 3 yr~', van den Bergh & Tammann ll991l Cappellaro & 
Turatto [19971 Li et al. l20TTb . For the case of a CE = 3.0, the 
SNe la start at the age of several tens of Myr and the birth rate 
increases remarkably until about 0.2 Gyr. The birth rate then 
slowly increases, again because SNe la with short delay times 
predominate the case of acE = 3.0. However, for the other two 
cases, the birthrates permanently increase with delay time, ow- 
ing to most of SNe la having a long delay time. 

Fig. [4] presents the effects of metallicity on the birth rate 
of SNe la with a constant SFR of 5M yr _1 . From the figure, 
we can see that decreasing the metallicity from 0.02 to 0.0001 
may increase the present birth rates by up to about 50%, which 
is qualitatively consistent with the discovery that the SN la rate 
in lower-metallicity galaxies is higher than that in metal-rich 
environments ( Kister et al. 1201 11 1. However, the present birth 
rate does not monotonically increase with decreasing metallic- 
ity, i.e. the highest value is for the case of Z = 0.002. The rea- 
son for this phenomenon is that explained in the above section, 
i.e. a low metallicity may lead to a longer GWR timescale. If 
the GWR timescale of DD system is so long that its delay time 
is longer than 15 Gyr, it does not contribute to the present birth 
rate of SNe la. Interestingly, the DD system was proposed to 
be the progenitor of super-luminous SNe la (Astier et al. 2006 ; 
Howell et al. 2006 ), which tend to explode in metal-poor envi- 
ronments (Khan et al. 12011b . 



the CE evolution, an earlier explosion is expected. Another fea- 
ture in Fig. |2]is that low metallicity generally leads to a SN la 
with a long delay time, especially for those older than ~ 1.5 
Gyr, which results in a shallower power law at low metallicity. 
This also reflects the effects of metallicity on the stellar evolu- 
tion. Generally, a star with a low metallicity usually produces a 
more massive white dwarf, i.e. a less massive CE (Meng, Chen 
& Han 2008). After CE ejection, the surviving DD system of 
low metallicity has a larger orbital separation than that of high 
metallicity. Although a more massive WD also has a shorter 
GWR timescale, the separation is the dominant factor when 
determining the GWR timescale, and a large orbital separation 
then means a longer GWR timescale. 

In Fig. |2] we compare our power-law DTDs with recent 
observations. Our results are marginally consistent with obser- 
vations. 



3.3. The normalization of DTD 

Apart from the form of the DTD, there is also fairly good agree- 
ment in terms of the DTD normalization Nsn/M*, i.e. the time- 
integrated number of SNe la per stellar mass formed. The range 
of its value is (0.5 - 3.5) x lO- 3 /M , focusing on 2 x 10" 3 /M o 
(Maoz & Mannucci l201 11 1. In Tab.[T] we present the DTD nor- 
malization for DD mergers with different metallicities and dif- 
ferent treatments of CE evolution. The results here are again 
marginally consistent with those from observations. The upper 
limit to the DTD normalization here is only 0.87 x 1O~ 3 /M , 
which is true for the case of Z = 0.002 and <?ce = 3.0. 

3.4. The DTD from SD + DD 

In Fig. |5j we show a DTD obtained from a combination of 
the SD and DD scenarios, where the SD scenario includes the 



Meng & Yang: The birth rate of supernovae from DD and CD systems 



7 



SD+DD - 




8 8.5 9 9.5 10 10.5 

log(t/yr) 

f ' SD+DD - 




J i l_j i i I i i_l i i L_ 

8 9 10 

log(«/yr) 



Fig. 5. The evolution of the birth rate of SNe la for a single star- 
burst of 1O U M and Z = 0.02. Solid, dashed, and dot-dashed 
lines are from SD + DD, SD, and DD channels, respectively. 
The results for the SD channel are from Wang et al. (2009bi 
and Meng & Yang d2010al > including the WD + He star, WD + 
MS, and WD + RG channels, while those for the DD channel 
are from this paper. The circle is from Mannucci et al. (12005) . 
the filled circles are from Totani et al. (12008) . stars are from 
Maoz & Badenes (2010), squares are from Maoz et al. (120 10) . 
and triangles are from Maoz et al. (1201 1) . Top: uqe = 1-0; 
Botom: acE = 3.0 

WD + He star, WD + MS, and WD + RG channels and the 
mass transfers for these channels all occur via a stable RLOF 
between a WD and helium/main sequence/redgiant star, not 
by wind accretion (see Wang et al. [2009b and Meng & Yang 
I2010al in details). For the WD + He star channel in Wang et al. 
(2009b), a helium star is in helium main sequence or crosses 
the helium HG. As in the DD channel, the DTD from the SD 
channel is also affected by the treatment of CE, so is the com- 
bined DTD. For oqe = 10, the SD channel dominates the SNe 
la with a short delay time, while the DD channel contributes 
most of the SNe la older than about 2 Gyr. The final combined 
DTD is a weak bimodality, where the WD + He star channel 
mainly contributes to SNe la with delay times shorter than 0. 1 
Gyr, WD + MS to those of 0.1 - 2 Gyr and DD channel to 
those older than 2 Gyr. At delay time shorter than 2 Gyr, this 
combined DTD is not inconsistent with observations for large 




20-156 







Fig. 6. The evolution of the distribution of the total mass of 
the DD systems leading to SNe la for a single starburst, where 
Z=0.02 and acE - 3.0. In the figure, all the cells move down- 
ward artificially by 1.378M . 



observational errors, but the long delay-time part of the DTD 
is significantly lower than observations by a factor up to 10. 
For ccE = 3.0, the combined DTD is closely consistent with 
the observational data for SNe la with delay times shorter than 
2 Gyr, while for the SNe la older than 2 Gyr, the combined 
DTD only marginally matches observations. As in the case of 
£*ce = 1 0, the DD channel is also the dominant channel for 
the SNe la older than 2 Gyr. The contribution of the WD + 
RG channel could be neglected (see also Lipunov et al. 1201 It . 
However, for young and middle age SNe la, the contributions 
of the SD and DD scenarios are comparable. 

3.5. The distribution of total mass 

Howell (1201 II) suggested that the DD scenario has a natural ex- 
planation for a relatively high SN luminosity in young environ- 
ments, namely that younger, more massive stars produce more 
massive white dwarfs, which have more potential fuel than less 
massive mergers, and then brighter SNe la (see also Maoz & 
Mannucci 1201 II ). If this suggestion were right, we should see a 
significant change in the total mass range of DD systems and 
its average value, i.e. the range and intermediate value of the 
total mass should decrease with delay time. In Fig. [6] we show 
the evolution of the distribution of the total mass, but we do 
not find any significant expected evolution, i.e. the distribution 
of the total mass is rather uniform across the whole delay-time 
interval. The uniform distribution is produced mainly because 
the dominant delay time for SNe la from DD systems is deter- 
mined by the GWR, not the evolutionary time of the secondary. 

We investigated various physical properties when search- 
ing for a correlation between delay time and a mass-dependent 
quality (the masses of more massive WD M\ , less massive one 
M2, mass ratio M2/M1, reduced mass u * +] £ , the radius of less 
massive WD, the final separation when coalescence begins, and 
the total angular momentum of the system at the moment of the 
onset of coalescence). Despite our efforts, no correlation fulfills 
the observational limit (the range of the parameter decreases 
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with the age of the SNe la, while the average value of the pa- 
rameter decreases/increases with the age) was found, which has 
roots in the evolution of the total mass and mass ratio distribu- 
tions. Hence, under the framework of DD scenario, it may still 
be difficult to explain the scatter in the brightness of the SNe la 
or the Phillips relation and more efforts are needed to build a 
bridge between the DD model and explosion model of SNe la. 

In Fig.|6l most of the DD systems have a mass higher than 
the Chandrasekhar mass by no more than 0.5 M Q . During the 
merging process of a DD system, the less massive WD is gen- 
erally destroyed and its mass is accreted onto the more massive 
WD. On the basis of Shen et al. (2012), we could expect that 
after the destruction of the less massive WD, a giant-like struc- 
ture be produced during the accretion stage. Because of the rel- 
atively long-lasting time that the giant-like structure persists, 
we might expect to lose about half a solar mass from the sys- 
tem (Soker l201 11 1. If this were a real scenario, according to the 
results in Fig. [6] the final masses of the most mergers could not 
reach the Chandrasekhar mass, and the SNe la could not then 
be expected, even though the total mass of the DD systems be- 
fore merging exceeds the Chandrasekhar mass limit. 

4. The birth rate of SNe la from CD systems 

4.1. Without spinning-down time 

In Fig.|7J we show the evolution of the birth rate of SNe la from 
the CD channel for a single starburst and a constant star forma- 
tion rate with different treatments of CE, where the spinning- 
down time from magneto-dipole radiation is assumed to be in- 
significant and SN la occurs immediately after merging. The 
cutting time means that only when a CD system merges within 
the cutting time, is it assumed to explode as a SNe la. In the 
figure, almost all of the SN la produced by CD systems have 
a delay time shorter than 2 x 10 8 yr and the birth rate peaks at 
~ 1.5 x 10 s yr if we assume that all the potential CD systems 
must merge within lOx 10 5 yr after the CD systems form, where 
the delay time is mainly determined by the evolutionary time 
of the primordial secondary whose primordial mass is gener- 
ally higher than 3 M . Hence, if a spinning-down time follow- 
ing magneto-dipole radiation did not play a dominant role on 
the delay time of SNe la, the CD scenario may only produce 
SNe la with very short delay times as expected by Livio & 
Riess ( 120031 1. In panel (1) of Fig.|7J some SNe la may have de- 
lay times as long as 1.5 Gyr, especially for ycE = 1-5, which 
is from the "pollution" channel. Our results for cutting times 
shorter than 10 x 10 5 yr are not polluted by the channel. 

Even in the case of f cut < 10 x 10 5 yr, the results are still 
significantly affected by the treatment of CE evolution. For a 
a-formalism, if acE < 2.0, no CD system survives the CE evo- 
lution. Even for oce ^ 2.0 (acE^ ^ 1-0), the CD systems that 
survived CE evolution will always merge within a timescale 
longer than 1 x 10 5 yr. In contrast, the y-algorithm may pro- 
duce CD systems with various merger timesc ales. 

However, the peak value of the birth rate for the CD channel 
is lower than that for the normal DD channel by 1 - 3 magni- 
tudes, which means that the SNe la produced in the CD sce- 
nario are not the main contributors to all SNe la. As shown 
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Fig. 7. The evolution of the birth rate of SNe la from a CD 
channel for a single starburst of 10 n M G (Top) and constant 
star formation rate of 5M yr~' (Bottom) with different cutting 
times and different treatments of the CE, where Z = 0.02. The 
delay time in the figure does not include the spinning-down 
time from magneto-dipole radiation. Solid, dashed, and dot- 
dashed lines are the cases of acE = 2.0, qtce = 3.0, and jce 
1 .5, respectively. 



in the bottom figure of Fig. [7] no more than 1 in 1000 SNe 
la should be produced by the CD scenario. The low birth rate 
results from the constraint that the more massive hot core orig- 
inates from the primordial secondary. 



4.2. With spinning-down time 

Ilkov & Soker ( 120121 1 suggested that the delay time for SNe la 
produced by a CD system may be dominated by the spin-down 
timescale related to magneto-dipole radiation torque, and the 
timescale from an initial fast rotation Qq to a critical angular 
velocity D. c is 



v -2 



> o.7n K 



^Uoookmy v o.i ; \o.3) 1 \n c ) 



(5) 



where B is the magnetic field, R is the radius of WD, and ji\ is a 
structural constant for the moment of inertia. For simplicity, we 
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Fig. 8. The evolution of the birth rate of SNe la from the CD 
channel for a single starburst of lO n M with different cutting 
times and different treatments of CE, where Z = 0.02. The de- 
lay time in the figure includes the spinning-down time from 
magneto-dipole radiation. Solid, dashed, and dot-dashed lines 
are the cases for acE = 2.0, oce = 3.0, and jqe - 1.5, respec- 
tively. 



assume that Q = Q Kep , Q. c - 0.7Q Kep , B\ - 0.3, and R = 4000 
km. The spin-down timescale is then 



t b ~ 5 x 10 7 



B sin 5 
10 7 G 



-2 



yr, 



(6) 



where B sin 5 follows a distribution of 



dN 



d\og{B sin 5) 



— constant, 



(7) 



for 10 6 G < B sin (5 < 10 8 G (Ilkov & Soker|20T2]). On the basis 
of this distribution, we obtained the spin-down timescale us- 
ing a Monte Carlo simulations, and inferred the delay time of a 
SN la to be the sum of the evolutionary timescale in forming a 
CD system, the one for GWR, and the spin-down timescale if 
the total mass of the CD system is lower than 1.48 M other- 
wise the spin-down timescale is neglected for those with total 
masses higher than 1.48 M (Yoon & Langer l2004l ). 

In Fig. [8j we show the DTD of SNe la produced by a CD 
system, where the spin-down timescale is included. In the fig- 
ure, the DTD also follows a power law of t~', where s ranges 
from 1.0 to 1.5 depending on the treatment of CE evolution. 
This is a natural result since the distribution of the B sin 6 is 
deduced from a power-law DTD (Ilkov & Soker 2012 1. Hence, 
the spin-down timescale is very important to the SNe la pro- 
duced by the CD channel since the power-law shape does not 
occur according to the a-formalism if the spin-down timescale 
is not taken into account (see Fig. [7]). In addition, as expected 
by Kashi & Soker d20TTb and Ilkov & Soker d20T2b . a few SNe 
la produced by CD systems may have delay times as long as 
several Gyr. 



5. Discussions and conclusions 

5.1. Other possible channels for SNe la 

From our study, we have found that the birth rate of SNe la 
for DD channels may only be marginally consistent with ob- 
servations, and the rate depends significantly on the treatments 
of CE. Only when oqe ^ 1-5, may we obtain a DTD whose 
shape follows a power law of t , and the birth rate is only 
consistent with the lower limit of this trend derived from ob- 
servations. For the cases of otce < 1-5 and ycE = 1-5, neither 
the DTD shapes nor the birth rates match the available obser- 
vations. Badenes & Maoz ( 120 121) calculated the merger rate of 
binary WDs in the Galactic disk based on observational data 
in the Sloan Digital Sky Survey and concluded that the merger 
rate of binary WDs with super-Chandrasekhar masses would 
not control significantly the SNe la rate. Our results are qual- 
itatively consistent with the calculations of Badenes & Maoz 
(1201 21) . If we consider the contribution of the SD channel as 
done in Mennekens et al. ([3010), the combination DTD may 
match observations for delay times t < 2.5 Gyr, at least within 
observational errors, regardless of the value of acE- However, 
when t > 2.5 Gyr, the birth rate from combination DTD is 
lower than that from observations. Therefore, other channels or 
mechanisms may contribute to SNe la with long delay times. 

In this paper, the WD + He star, the WD + MS, and the 
WD + RG channels are included for the SD channel, while the 
wind-accretion channel is exclued. Chen, Han & Tout (1201 11 ) 
constructed a tidally-enhanced wind-accretion model where the 
initial mass of RG donor may be as low as 1 M , i.e. the delay 
time from the channel may be as long as 10 Gyr. On the basis of 
Eq. (1) of Iben & Tutukov ( 1984), they obtained a birth rate of 
6.9 x 10- 3 yr~' fr om the wind model. Although we believe that 
the birth rate is probably overestimated since some parameter 
spaces considered to be valid for SNe la production in Eq. (1) 
of Iben & Tutukov d 19841 ) may not contribute to SNe la, the 
wind channel should improve our results in this paper. 

Some overluminous SNe la were observed and a super- 
Chandrasekhar mass explosion was expected for these SNe la 
based on the amount of 56 Ni inferred (Astier et al. 2006 Howell 
et al. 120061 Hicken et al. 120071 Scalzo et al. 120101 Yuan et al. 
l20T0l Tanaka et al. l20T0l Yamanaka et al. l20T0l ). These overlu- 
minous SNe la are most likely to occur in metal-poor environ- 
ments (Khan et al. 1201 lb . Although they are generally assumed 
to be produced by the coalescence of DD systems, they may 
also originate in SD systems (Chen & Li 120091 Liu et al.|2009; 
Hachisu et al. 2012 ). However, the super-Chandrasekhar mass 
explosion from the SD channel may only contribute to SNe la 
with short delay times (Hachisu et al. 2012), and their contri- 
bution to all SNe la may be no more than 0.3% (Meng et al. 
1201 m . 



5.2. The effect of metallicity 

As noted in Sects. l3~D and !3. 21 the birth rate of SNe la depends 
on metallicity, especially for those with long delay times. Meng 
et al. ( 1201 lal ) found that decreasing the metallicity may signif- 
icantly increase the birth rate of SNe la with long delay times, 
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that are produced by the WD + RG channel, by a factor of about 
three. Some observations do indeed find that a lower metallic - 
ity may correspond to a higher SNe la birth rate (Kistler et al. 

The picture emerging from some observations is remark- 
able, i.e. most of these diverse DTD derived from different 
methods, different environments, and different redshifts agree 
with each other, in both form and absolute value (Maoz & 
Mannucci fcOllI Graur et al. 1201 11 1. At delays of t > 1 Gyr, 
a power law of index about - 1 seems to be in little doubt, al- 
though a index of ~ -1.5 could not be ruled out (Barbary et al. 
120121 Sand et al. 120121 ). while at delays t < 1 Gyr, the shape 
of the DTD might become either shallower or steeper, or fol- 
low the same shape seen at long delays. In this paper, we note 
that for the DD scenario, a low metallicity may lead to a shal- 
lower DTD. Interestingly, for a SD channel, a shallower DTD 
may also be obtained at low metallicity (Meng et al. 1201 lab . In 
addition, metallicity may also affect the value of the birth rate 
of SNe la (Kistler et al. 1201 lb . Nevertheless, the level of the 
effect of metallicity on the DTD might not be as high as we ex- 
pect (see in Fig. [2]) since the DTDs from different environments 
agree with each other. However, we still have to determine us- 
ing observations the precise effect of metallicity on the DTD 
since metallicity does have an influence on the DTD, on both 
its form and absolute value, and then the shape of the DTD at 
delays t < 1 Gyr might become clear. 

5.3. The origin of the Phillips relation 

The Phillips relation is the most fundamental relation when 
a SN la is used as a distance indicator, which implies that 
the properties of SNe la are mainly dominated by one param- 
eter. However, the nature of this parameter remains unclear 
(Podsiadlowski 2008 ). Many efforts have been made to resolve 
this problem, but most of these discussions have focused on the 
Chandrasekhar mass model, in which the WDs explode as SNe 
la when their masses are close to the Chandrasekhar mass limit. 
No consensus has yet been reached at present. Here, we tried to 
identify a parameter in the DD scenario to explain the scatter in 
the luminosity of SNe la. Unluckily, no such a parameter was 
found. In particular, the total mass of the DD systems is a poor 
choice of parameter to illustrate the dependence of the average 
luminosity of SNe la on their age. Hence, the origin of Phillips 
relation still remains unclear. 

5.4. The birth rate of SNe la from CD channel 

We have found that the upper limit to the contribution of the 
SNe la produced by the CD channel to the total SNe la is 
only about 0.1%. One may argue that this result may seriously 
depend on the treatment of CE evolution since the treatment 
here is very simple and the real scenario may be very complex 
(Kashi & Soker l20TT] Passy et al. 2012). At present, it remains 
very difficult to construct a more complete model in a BPS 
study. However, since the simple treatment may succeed in re- 
producing the forms of many special objects, such as planetary 



nebulae and subdwarf B stars (Han et al. [T9951 120021 120071 ). 
the results here could be reasonable to at least first order. 



However, one may remark that since the CD systems had 
hardly survived the CE evolution, we might ignore the systems 
that almost survive the CE evolution, especially for a lower 
(ten- These systems might also explode if their spinning-down 
time is long enough to permit the remaining CE material to be 
lost by a wind. Because the results of 0.1% is derived from 
the a-formalism of oce ^ 2.0, we might underestimate the 
birth rate of SNe la predicted by the CD scenario. However, 
if q-ce < 2.0, no CD system can survive the CE evolution. 
This is a natural result since a low otce means that a system 
needs to shrink its separation more dramatically and release 
more orbital energy to eject the CE. We do not verify whether 
this merger occurs because this would require consideration of 
the details of CE evolution and many assumptions have to be 
made, which may lead to too many uncertainties. Fortunately, 
our results place a constraint on the problem. For a high o;ce, 
the systems that merge under the situation of low o;ce could 
survive, while those surviving the CE evolution with low otce 
may lead to a larger orbital separation. Hence, the panel (1) in 
Figs. Q and [8] should give a safe upper limit on the birth rate 
of SNe la based on the CD scenario, where all the potential 
CD systems should be included. Even so, the total contribution 
from all the potential CD systems to all SNe la is no more than 
1%. In addition, the spin-down mechanism is very important 
to the DTD shape of SNe la derived for the CD systems, 
especially for the a-formalism. We cannot obtain a power-law 
form if the spin-down timescale due to the magneto-dipole 
radiation torque is not taken into account. 

In summary, we have calculated the evolution of the birth 
rate of SNe la from DD and CD channels. We found that the 
treatment of the CE evolution has a great influence on the final 
DTD shape, and only for the a-formalism when oce ^ 1-5, 
the shape of the DTD is consistent with that derived from 
observations, but the birth rate of this case only marginally 
matches that derived observationally. For a cr-formalism with 
a lower c?ce and the y-algorithm, the birth rate of SNe la 
from DD systems is much lower than those derived from 
observations, and the shape of the predicted DTD does not 
follow a power law. Metallicity has almost no influence on 
the shape of DTD except in terms of the index of the power 
law, but may increase the value of the birth rate of SNe la. 
When the SD channel including WD + He star, WD + MS, 
and WD + RG channels is incorporated, the theoretical DTD 
matches the observations very well for SNe la younger than 
2.5 Gyr, while for those older than 2.5 Gyr, the theoretical 
birth rate is slightly lower than that derived from observations. 
We also need to incorporate other channels or mechanisms 
that contribute to the SNe la explosion. As suggested by 
Soker (120111 ) and Ilkov & Soker (120121 ). the CD scenario is 
a possible channel contributing to SNe la and we have found 
three potential channels that may produce the CD systems. 
However, the strict upper limit to the contribution from the CD 
scenario to all SNe la is 1%. 
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